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The new data on the 7*7— >tt° transition form factor of the Belle Col- 
laboration are analyzed in comparison with those of BaBar (including the 
older data of CELLO and CLEO) using an approach based on light-cone 
sum rules. Performing a 2-, and a 3-parametric fit to these data, we found 
that the Belle and the BaBar data have no overlap at the la level. While 
the Belle data agree with our predictions, the Babar data are in conflict 
with them. 

PACS numbers: 12.38.Lg, 12.38.Bx, 13.40. Gp, ll.10.Hi 

1. Light-cone sum rules for the process 7*(Q 2 )7(q 2 — 0) — > tt° 

The validity of the collinear factorization — the basis of applications of 
QCD to hard processes — was challenged in the year 2009 by the exper- 
imental data measured by the BaBar Collaboration pQ for the kinematics 
Q 2 > m 2 q 2 <C m 2 . This year, new experimental data by the Belle Collab- 
oration [2] were presented that do not indicate such a growth and are grossly 
compatible with the QCD expectations. In this presentation, we show how 
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these data and the previous ones by the CELLO |3] and CLEO [1] Collab- 
orations can be analyzed within the theoretical scheme of light-cone sum 
rules (LCSR)s that incorporates contributions from QCD perturbation the- 
ory and higher-twist corrections. Within QCD, the pion-photon transition 
form factor (TFF) 7 71 is given by the matrix element 



d 4 ze~^ z (7r°(P)\T{j ll (z)j u (0)}\0) 



H 



r, a r,P T?~l*~l*-K° m 2 r, 2 



(1) 



where denotes the quark electromagnetic current and Q 2 = —q\ > 0, q 2 = 
—q\ > 0. For the asymmetric kinematics q 2 < m 2 , one has to include into 
the calculation the interaction of the (quasi) real photon at long distances 
~ 0(1/ y/q 2 ). To accomplish this goal, we apply the approach of LCSRs 
[6] that effectively accounts for the effects of the long-distance interactions 
of the real photon by making use of a dispersion relation in q 2 and applying 
quark-hadron duality in the vector channel [5]. Taking the limit q 2 — > 0, 
we get [B] 



Q 2 F^° (Q 2 ) = ^ 



Q 2 



771, 



exp 



/ P(Q\x)- 



Q 2 x/x 



M 2 



p(Q 2 ,x) 



dx 



(2) 



where M 2 is the Borel parameter in the interval 0.7-0.9 GeV 2 and the spec- 
tral density is given by p(Q 2 ,x) = (Q 2 + s)p pert (Q 2 , s) with 



P 



pert 



1 



(Q 2 , s) = -ImF 7 * 7 * 71 " (Q 2 , -s - ie) = p t w-2 + Ptw-4 + Ptw-6 + 



(3) 



The various twist contributions are defined in the form of convolutions of 
the corresponding hard parts with the pion distribution amplitude (DA) of 
a given twist [6]. For instance, for the twist-four contribution we use the 

effective description [6] ip^\x,fi 2 ) ~ <5 2 w _ 4 (/x 2 ) x 2 (l - x) 2 with (5 2 w _ 4 (/j 2 ) = 
0.19 ±0.04 GeV 2 [7J. Here we used the abbreviations x = 1 — x, s = xQ 2 /x, 
%o = Q 2 /(Q 2 + sq), where sq ~ 1.5 GeV 2 is the effective threshold in the 
vector channel. The leading twist-two contribution has the perturbative 
expansion (a s = a s /(47r)) 



ritw-2 

-T * » n ~ 

7*7*7T U 



(T LO + a s (// 2 )T NLO + a 2 (/i 2 )T NNLO/3 + . . .) ® ^(x, fj 2 ) . (4) 



The corresponding contributions from Q to the spectral density §3§ 
have been obtained in [§]. For the term pnlo we employ the corrected 
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version computed in [lQj . The "bunch" of the admissible twist-two pion 
DAs </? T (x,// 2 ) was determined in [11] within the framework of QCD SRs 
with nonlocal condensates (NLC)s. This DA "bunch" (shown graphically in 
Fig. [1] in the form of a rectangle) can be effectively parameterized in terms 
of the first two Gegenbauer coefficients ai and (24 to reduce the expansion to 



Our next goal is to 



<p™ s (x) = Qx3 , 1 + a2 c^ 2 (2x - 1) + a 4 Cl /2 {2x - 1) 

extract a pion DA that provides best agreement with all sets of the existing 
data by performing a fit procedure of the Gegenbauer coefficients a n within 
the basis of LCSRs, cf. Eq. ([2]). The results for the combined set of data 
from CELLO&CLEO&Belle - termed CCBe — will be contrasted to those 
from the CELLO&CLEO&BaBar set — called CCBB. 



2. 2D-analysis of the combined data fit 

In Fig. [TJ we present the confidence region of the 02, 04 values in the form of 
error ellipses in the (02,04) plane, obtained by fitting different sets of data. 
We take into account only statistical errors and exclude theoretical uncer- 
tainties — in contrast to our previous work in [TJ Ej- Using LCSRs with 
the pion DAs obtained in QCD SR NLC [11], one arrives at the predictions 
shown in terms of the (green) rectangle in comparison with the error ellipses 
pertaining to the different sets of data defined in the previous section. The 
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Fig. 1. Error ellipses of various sets of data described in Sec.[T]and taken at fj, = 2.4 GeV 
scale. The best-fit values of \ 2 are shown at the centers of the ellipses, while differences 
between data sets are displayed along the links in units of la. The vertical lines show the 
range of values computed on the lattice: dashed line — [12]; dashed-dotted line — [13] . 

best-fit values of the x 2 goodness criterion xidi = X 2 / n df (ndf= number of 
degrees of freedom) are shown as centers of the ellipses, with the devia- 
tions from one data set to another with reference to the rectangle being 
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displayed along the links and expressed in units of one standard deviation 
(la ~ 68%). We conclude that the inclusion of the BaBar data to those 
obtained before by CELLO&CLEO leads to an approximately 3a shift [13] 
of the confidence region away from the (black) ellipse at the bottom. The 
result is represented by the (red) dotted ellipse at the top accompanied by 
a significant increase of xidi f rom 0-4 to 2 (for the CCBB data set). If we 
include to the CELLO&CLEO data only the Belle data set CCBe, then the 
shift of the confidence region is only moderate giving rise to a slight increase 
of xidf f rom 0-4 to 0.6 (ellipse in the middle) \j We quantify the deviations 
of these data sets from our theoretical estimates, encoded in the rectangle, 
in terms of a values shown along the links which connect the central data 
values with the BMS model (X) inside the rectangle. Here, and in the next 
figure, the vertical broken lines denote the constraints from two different 
lattice simulations: dashed lines — [12]; dashed-dotted lines — [13j . 
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Fig. 2. (color online). Predictions for the moments (£ 2 ) and (£ 4 )tt at the lattice scale 
MLat — 4 GeV 2 . The lower stretched ellipse — solid (blue) line — corresponds to a 2D fit 
to the CCBe data, while the upper one (red line) represents the CCBB one. In contrast 
to Fig. [TJ we also include here the uncertainty related to the twist-four coefficient (5 2 w _4. 
The vertical lines show the range of values computed in two lattice calculations using the 
same designations as in Fig. [T] 

An immediate conclusion from these considerations is that (i) the BMS 
DA is inside the la CCBe and inside the 0.6a region of the CELLO&CLEO 
set, whereas the BMS "bunch" greatly overlaps with both of these error 
ellipses. At the same time, the central point of CCBB is 6.2<j away and 
has no overlap with the CCBe ellipse, (ii) The asymptotic DA (♦) and the 
Chernyak-Zhitnitsky (CZ) DA (■) are both more than 6a away from CCBe. 
(iii) The existing lattice calculations of 02, shown as vertical lines in Figs. 

1 Note that we consider the x 2 goodness criterion for the considered experiments as 
being the sum of the individual % 2 values associated with each experiment. 
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11121 are not restrictive enough to narrow down the interval of 02, though 
the narrower band [13] supports both the CCBe data and our theoretical 
results but not the CCBB data. On the other hand, the previous lattice 
simulation [12] is compatible with all sets of data and with our theoretical 
predictions as well. 

With all these data in our hands, we may attempt to extract values of 
the moment (£ 4 )tt by combining the real data from CCBe and CCBB with 
the results of lattice simulations. The combined constraints from CCBe 
(lower (blue) error ellipse in Fig. [2]) and the lattice constraints from [13] 
(dashed-dotted vertical lines) lead to the following predictions for (£ ^ 
li 2 )^ £ [0.26 0.30] and (i A ) n G [0.112 0.13]. In contrast, the constraints 
extracted from the CCBB data (upper (red) ellipse) do not allow a reliable 
determination of the (£ 4 )tt range in the background of the lattice results 
from |13j . On the other hand, the constraints from [12] are not stringent 
enough to differentiate these sets of data. 

3. 3D-analysis of the combined data 

We step up from a 2D to a 3D analysis of the data, by including the 
next coefficient ciq. Then, we get in Fig. [3] fit results in terms of la ellipsoids 
with respect to the experimental statistical errors for the CCBe data set 
(left panel) and the CCBB (right panel). The theoretical =pA(5^ w 4 -error is 
indicated by a solid (red) cross for an increasing value of of w _ 4 , whereas a 
dashed (green) cross (closer to the a& axis) denotes a decreasing value. The 
projection of the CCBe la ellipsoid on the plane (02,04) is represented, in 
both panels, by the larger (red) ellipse, while the smaller one refers to CCBB. 
The shaded (green) rectangle encloses the region of 02,04 pairs allowed by 
NLC SRs p], with the symbol (*) marking the BMS pion DA. All results 
are shown at the scale ji = 2.4 GeV after NLO evolution. 

To further quantify these statements, we supply the "coordinates" of the 
central point of each of the two displayed ellipsoids in the following form. 
The first number gives the central value of the fit, the next number is the 
statistical la error, and the third one is the theoretical uncertainty due to 
twist-four. Then, we have: CCBe (0.157 ± 0.057 ± 0.056, -0.192 ± 0.122 ± 
0.077, 0.226 ± 0.161 ± 0.033) with x„df = 13.1/30; CCBB (0.177 ± 0.054 ± 
0.056, -0.171 ± 0.103 ± 0.071, 0.307 ± 0.096 ± 0.024) with Xndf = 33.3/32. 
To conclude: (i) The description of the CCBe data provides a better x„df ~ 
0.4 relative to x„ d f ^ 1 following from CCBB. (ii) The CCBe and CCBB 
ellipsoids are significantly separated from each other, so that a QCD de- 
scription of these sets of data requires substantially different DAs. (iii) The 
(02,04) projections of both ellipsoids have a good overlap with the BMS 
"bunch", though the CCBe ellipsoid that has no intersection with it. An 
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Fig. 3. 3D graphics of the la error ellipsoids in the space spanned by the Gegenbauer 
coefficients a n (n = 2, 4, 6) for processing the data on the pion-photon transition form 
factor within the LCSR approach using the scale /x=2.4 GeV. Left. Only CCBe data 
used. Right. CCBB fit (smaller ellipsoid) shown in comparison with CCBe data. More 
explanations are given in the text. 



intersection is possible but only at a larger value of confidence level. 



4. Local characteristics of the pion DA 

The confidence region of the coefficients {a n }, obtained in Sec. can 
be linked to other characteristics of the pion DA. The profiles of the pion 
DA (p v (x), extracted in the 3D fit procedure, are shown in Fig. left panel 
- set CCBe; right panel — set CCBB. The BMS "bunch" (shaded green 
strip) and the BMS DA model (black solid line inside it) are also shown. 
The main difference between the two graphics in Fig. [5] stems from their dis- 
tinct behavior near the endpoints that are concentrated inside the interval 
A = 0.05. Indeed, one observes that the DA (including uncertainties) - 
shaded area in the right panel — providing best-fit to the CCBB data devi- 
ates significantly from the BMS "bunch". A mathematical tool to quantify 
the endpoint characteristics of the pion DA is L>( m V(A), i.e., the average 
derivative of ip(x) in the interval A, that was invented in [TS] and possesses 

the following properties: lim D^ip n (0)(A) = (p'(0); lim Z)( m )^(0)(A) = 

A— K) m— ><x> 
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Fig. 4. Left. Comparison of the BMS pion DA "bunch" (darker shaded strip in green) and 
of the BMS model (black solid line inside this strip) with the 3D fit to the experimental 
data on the pion-photon transition form factor. The solid blue line denotes the best-fit 
pion DA (at the scale fi = 2.4 GeV, see second paper in [6]) obtained from the analysis of 
the set CCBe, with the dashed lines indicating the sum of the statistical errors of the fit 
and the twist-four uncertainties. Right. Analogous results obtained with the set CCBB. 



ip' w (0), while D^ip(l) = dx(p n (x)/x. The quantity D^ip(A) can be used 
[16| to probe the endpoint behavior of these different sorts of DAs: 

£> (2) y?(0.05) = 17.2 ± 8.5 for CCBe; 25.6 ± 5.25 for CCBB . (5) 

We see that the slope in the endpoint region of the best-fit DA to CCBB 
is much stronger than for CCBe. While the CCBe-based DA profile agrees 
with the BMS "bunch", it is incompatible with the CCBB shape. 



5. Conclusions 

We presented here an analysis of all experimental data on the pion- 
photon transition form factor and compared in detail the results obtained 
with two different sets of data: CELLO, CLEO, BaBar — (CCBB) vs. 
CELLO, CLEO, and Belle — (CCBe). Our analysis is based on LCSRs at 
NLO, also taking into account the twist-four term. The NNLO^ radiative 
correction was included into the theoretical uncertainties together with the 
twist-six contribution [THj- The key results can be summarized as follows: 
(i) We performed a 2D (parameters 02,04) and a 3D (parameters 02,04,06) 
analysis, fitting both the CCBe and the CCBB data sets. We found that 
in both cases the fit to CCBe provides a significantly better x^df va l ue - (h) 
The CCBe 2D fit agrees well within error bars with the previous findings 
from the CELLO and CLEO fits and the predictions derived from QCD SR 
NLC[TT] with the BMS "bunch" and the BMS model DA. The fits are also 
compatible with the constraints on 02 extracted from lattice computations. 
In contrast, the CCBB 2D fit has no overlap with the CCBe result and it 
also does not comply with the QCD SR NLC predictions, (iii) The results 
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of the 3D fit to the CCBe and CCBB data sets do not intersect at the level 
of a la accuracy. At the same time the CCBe result is much closer to the 
BMS "bunch" , but still outside the la area, (iv) The qualitative features of 
the 3D fit to the CCBe and CCBB data can be differentiated in terms of the 
behavior of the corresponding profiles of the pion DAs near the endpoints: 
the CCBe has a less pronounced slope and is close to the BMS "bunch", 
whereas the slope of CCBB is larger and gives no support to the BMS DAs. 
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